Scaling of Low-Order Structure Functions in Homogeneous Turbulence 
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High-resolution direct numerical simulation data for three-dimensional Navier-Stokes turbulence 
in a periodic box are used to study the scaling behavior of low-order velocity structure functions 
with positive and negative powers. Similar to high-order statistics, the low-order relative scaling 
exponents exhibit unambiguous departures from the Kolmogorov 1941 theory and agree well with 
existing multiscaling models. No transition from normal scaling to anomalous scaling is observed. 
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The intermittency effect of small scale dynamics on 
fluid turbulence has been an important issue in turbu- 
lence statistics for many decades. A lot of effort has 
been invested towards the understanding of the depar- 
ture from Gaussian statistics and the related anomalous 
scaling behavior of velocity structure functions. At very 
high Reynolds numbers, Kolmogorov's 1941 similarity 
theory (K41) predicts a normal (or regular) scaling for 
moments of the velocity increment in the inertial range: 
{\l\Ur\'^) ~ r'^/^ for L » r >> 77, where the longitudinal 
velocity increment is defined as Am^. = u{x + r) — u{x), 
(•) denotes an ensemble average and L and 77 are, re- 
spectively, the integral scale and the Kolmogorov dissi- 
pation scale. For convenience, we have considered here 
the moments of the absolute values of velocity differences 
— called the "generalized structure functions" to distin- 
guish them from the classical ones. The normal scaling 
has been challenged often and the departure from K41 
has been reported from both experiments and direct nu- 
merical simulations A recent paper used data 
from a high-resolution simulation and demonstrated the 
existence of intermittency correction to K41 without in- 
voking the Taylor's hypothesis used in real-life experi- 
ments. 

It has been believed traditionally that the deviations 
from K41 tend to be more pronounced when the mo- 
ment order q becomes large. This is because the non- 
Gaussianity is more conspicuous at the tails of the prob- 
ability density function (PDF) which make increasingly 
significant contributions to higher order statistics. Thus, 
experiments and simulations have focused on structure 
functions with > 2, say. With the exception of Refer- 
ences § and 0, little attention has been paid to the 
scaling of low order structure functions {q < 2, say). 
The preliminary experimental results of a water flow in a 
pipe in Ref. Q showed that exponents for moment orders 
as low as 0.25 departed from K41. However, a detailed 
study connecting large amplitude events with structure 
functions and scaling exponents has not yet been made. 

For one-dimensional Burgers turbulence in the small 
viscosity limit, a bifractal structure of velocity due to 
the formation of shock structures leads to a normal scal- 



ing for g < 1 and a constant scaling exponent for g > 1 
Q . This bifractal scaling has also been reported for tur- 
bulent convection with a constant temperature gradient 
One cannot in principle rule out a somewhat simi- 
lar situation in three-dimensional turbulence at very high 
Reynolds numbers. It is therefore both interesting and 
useful to examine the nature of low-order scaling expo- 
nents. 

In this Letter, we present an analysis using data 
from direct numerical simulations for a fully developed 
isotropic turbulence to study low order moments of ve- 
locity structure functions and the relation between struc- 
ture function exponents and the PDF. The particular 
questions we address here are: (1) Are there intermit- 
tency corrections to the scaling of low order structure 
functions? (2) Is there a transition from regular scaling 
at low order to anomalous scaling at high order? (3) 
What aspects of turbulence make major contribution to 
low order structure function statistics? 

Direct numerical simulation of the Navier-Stokes equa- 
tions [0 was carried out with 512^ mesh points in a 
cyclic cubic box for homogeneous isotropic turbulence 
using the CM-5 machine at Los Alamos and the SP ma- 
chines at IBM. The simulation domain was [0, 2tt] in 
each direction. A nominal steady state was maintained 
by a forcing confined to wavenumbers /c < 3. The Tay- 
lor microscale A = (ISi^Wq/e)^/^ and microscale Reynolds 
number TZ\ = vqX/v were controlled by varying the vis- 
cosity of the system. Here Vq is the root-mean-square 
value of a single vector component of velocity. The anal- 
ysis was carried out for forced statistically steady states 
at TZ\ — 218. A spatial averaging over the whole physical 
space was used to replace the ensemble average. The sep- 
aration r was taken along the x direction and the PDFs 
and statistics shown below were averaged over 1.5 large- 
eddy turnover times. 

In the inset of Fig. 1, we present the generalized struc- 
ture functions 

S,{r) = {\Aurn^r'^^, (1) 

with q = 0.4 and 0.8. Here is the scaling exponent. 
A power law range can be identified for 0.2 < r < 0.8. 
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In Ref. [|6|, we have shown that the third-order structure 
function in this region has the expected inertial range be- 
havior [Q. The structure functions are compared with 
a phenomenological multiscahng model proposed by She 
and Leveque (sohd hues). To better calculate the 
scaling exponents, in Fig. 1 we have utilized the so-called 
extended scaling similarity (ESS) hypothesis advocated 
by Benzi et al. i.e., plotting Sq{r) against 6*3 (r), and 
identified the relative scaling exponent, Cq- From our pre- 
vious experiences Q and the current study, we believe 
that the ESS works well for velocity structure functions 
in homogeneous turbulence — meaning that the relative 
scaling shows a wider region than the direct scaling re- 
gion shown in the inset. We have studied many instances 
of direct and relative scalings for different power indexes 
I psf , ranging from q — —0.8 to 12, and found the scaling 
relations to be very similar to those presented in Fig. 1. 
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FIG. 1. Velocity structure function Sq{r) for q — 0.4 and 0.8 
versus S'3(r). The inset shows the velocity structure functions 
versus r. The symbols are direct numerical simulation results, 
and solid lines indicate slopes from the SL model. 

In Fig. 2 we have compared the ESS scaling expo- 
nents for —0.8 < g < 2 with predictions from K41, a log- 
normal model by Kolmogorov (K62) Q which gives = 
g/3-|-g(3— g)^/18 where ^ = 2/9 is the dissipation scahng 
exponent pql, and a log-Poisson model by She and Lev- 
eque (SL) which gives Cg = (7/9-h2/3[l-(2/3)«/3]. To 
be more quantitative, we list here five scaling exponents 
from direct numerical simulations: C-o.s = —0.32 ±0.01, 
C0.2 = 0.074 ± 0.002, C0.4 = 0.150 ± 0.002, C0.6 = 
0.223 ± 0.002, C0.8 = 0.296 ± 0.002. For very smaU q, 
C^q should be a linear function of q following a Taylor 
expansion. The difference of the slope of scaling expo- 
nents when q ^ between K41 and SL model (or K62) 
is about 12% It is evident that the scaling exponents 
deviate from K41 even for very small q, in favor of the 
intermittency models j^jl^,|l^ for both positive and neg- 
ative powers. Since the departure from K41 starts from 
very small q and the relative scaling exponents agree well 
with some intermittency models for q from —0.8 up to 8 



, there is no evidence to support a transition of scaling 
exponents and the bifractal structure, in contrast to the 
case of Burgers equation ||^ or pressure-less turbulence 
0- 
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FIG. 2. Scaling exponents of structure functions as a function 
of the order index q compared with various models. 

As pointed out in it is generally believed that the 
PDF of the velocity increment has a Gaussian core for 
small amplitude events and a stretched exponential shape 
for the tail part. The intermittency effect presented in 
Fig. 2 for low order moments may imply two possibilities: 
(1) the low order moments have substantial contributions 
from the high amplitude (tail) events; (2) the core part of 
the PDFs is also sufficiently non-Gaussian. To shed light 
on these possibilities, we have enlarged in Fig. 3 the core 
part PDF of Aw^ with the inset showing the whole PDF. 
While the deviation from the Gaussian is evident at the 
tails, even the core part of the PDF displays departures 
from Gaussianity. This departure may have significant 
contributions to low order statistics. 





1 ' 

:M 






1 ^ 




1 1 1 1 : 
















:/ 10"' 


■ 












/ 10"^ 


■ 












/ 10' 














/ 10" 














10 = 


■ 












10"' 














10' 














10"° 












' ' ' 


1 . 




-5 



1 


5 


, 1 , \'- , 



-2 2 

A,u/S,(r)^" 

FIG. 3. The core part of the PDF for velocity increment (solid 
line) compared with Gaussian distribution (dotted line) . The 
inset shows the whole distribution for the PDFs. All curves 
are normalized by the root-mean-square. 
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To probe this further, we define a cumulative structure 
function, 

Sg{r,<f>)= r \<f>rP{<l>')dcf>', (2) 
J-(t> 

and the cumulation ratio 

where (j)' = Aur/Srms{^) Srms{r) = 'S'2(r)^/^ is the 
root-mean-square for a given separation r; P{<f)') is the 
probability density function of 4>'; (j) is the bound of trun- 
cated (j)' domain. The ratio 7g(0) quantitatively defines a 
cumulative contribution from events with amplitude up 
to (j) for the g'^'-order structure function. In Fig. 4, wc 
show 7g(0) as a function of </> for a separation distance 
in the inertial range (r = 0.29). For fixed 7g (</>), we see 
that the larger the q the larger the (f), consistent with the 
fact that the peak in (ff^P{<l)) shifts to larger values of (j) 
for increasingly high-order moments. Wc provide in the 
inset the complementary data on (p (the inverse function 
of 7 in Eq. (3)) as a function of q when 7q((/>) = 0.9. 
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FIG. 4. The cumulation ratio 7q((^) as functions of (/> for dif- 
ferent q (from left to riglit): q = -0.8, -0.5, -0.1, 0.1, 0.5, 
1, 2, 4, 6, 8. The separation distance r = 0.29 lies within 
the inertial range. The inset shows values corresponding to 
the intersections of the curves by the horizontal line given by 
7<, = 0.9. 

The aspect just discussed has been made more specific 
in Fig. 5 for the flatness and skewness of the velocity in- 
crement (r = 0.29). The data are compared with those of 
a Gaussian field. The departure from Gaussian statistics 
for the flatness starts from 2± 1, which might there- 
fore usefully define the boundary between the core and 
the tail of the PDF. Both the skewness and the flatness 
converge to constants (3.6 and -0.32 respectively) when 
~ 4, consistent with the inset in Fig. 4. It is interesting 
to note that the skewness has a crossover from positive to 



negative at (/) = 2, indicating that the negative skewness 
basically comes from events of intermediate amplitudes 
between about 2 and 4 standard deviations. 




-0.4 ^ 1 

2 4 6 8 



FIG. 5. The skewness and flatness factor as function of for 
r = 0.29 calculated from the data (bullets), compared with 
Gaussian values (solid lines). 



In order to sec the effects of the PDF on scaling expo- 
nents, in particular on the variation of scaling exponents 
with the variation of the upper bound, (p, it is useful to 
rewrite Eq. (1) as Cq{r) = din Sq{r) /dlnr. We flrst note 
that the scaling exponents represent a relative variation 
of the structure function as a function of the separation 
distance. As in Eq. (2), we can now define the g'^'-order 
cumulative scaling exponent as 

dlnr 



In Fig. 6, we show Cqi'P'), the average of (q{r,(j)) over 
inertial range, as a function of for g = 0.4, 2 and 8. 
Two features in this plot are worth noting: First, for 
(/) ~ 1, the scaling exponents for all three cases are close 
to K41. Second, all three cases converge to the SL mul- 
tiscaling model when </> ~ 6. This result implies that — as 
far as the scaling exponents are concerned — there is no 
fundamental difference in the relation between structure 
function and the PDF of the velocity increment for high 
order and low order statistics. Note that the transition in 
Fig. 6 from K41 to the SL (or other) multiscale results is 
smooth, which seems to support the multifractal picture 
in which the correction to anomalous scaling comes from 
all amplitudes with different fractal dimensions. 
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FIG. 6. Scaling exponents (q{<j)) versus (j) for q — 0.4,2 and 
8 for truncated structure functions. The values for complete 
structure functions from K41 and SL are also plotted for com- 
parison. 

It is not often appreciated that the non-Gaussianity of 
PDF of velocity increments is related to the phase coher- 
ence among velocity mode in Fourier space ("structure" 
in real space): when a pseudo- velocity field is generated 
by randomizing the phase for each mode (but not altering 
its amplitude) the resulting scaling exponents agree with 
K41. Although this result is not surprising, it confirms 
the direct link of intermittency and anomalous scaling to 
real space structures 

To summarize, we have reported a study of low or- 
der scaling for the direct numerical simulation data on 
isotropic turbulence. The scaling exponents for low order 
structure functions clearly deviate from K41 but agree 
well with existing multiscaling models. No transition 
from regular scaling to anomalous scaling can be iden- 
tified. The dependence of structure functions and scal- 
ing exponents on amplitudes of the velocity increment 
are investigated. While the low and high order struc- 
ture functions have primary contributions from small 
and large amplitude events, respectively, the scaling ex- 
ponents for structure functions, regardless of the power 
index q, seems to have some contribution from all events. 
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